RECENT PROGRESS IN THE CALDERON PROBLEM 
WITH PARTIAL DATA 



CARLOS KENIG AND MIKKO SALO 

Abstract. We survey recent results on Calderon's inverse prob- 
lem with partial data, focusing on three and higher dimensions. 
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1. Introduction 

The Calderon problem with partial data asks to determine the elec- 
trical conductivity of a body from electrical measurements made on 
certain subsets of the boundary. This is a fundamental inverse prob- 
lem, also mentioned as an open question in Gunther Uhlmann's ICM 
address |Uh98] . Subsequent years have seen several advances in partial 
data problems, many of them due to Gunther Uhlmann and his collab- 
orators, and it is a pleasure for us to survey some of these developments 
in this volume in honor of Gunther's remarkable career. 

Let us give the formal statement of the problem. If fi C M™, n > 2, 
is a bounded domain with C°° boundary, if 7 e L°°(Q) is a positive 
function (the electrical conductivity of the medium), and if Tp, TV are 
open subsets of dQ, we consider the partial Cauchy data set 

C^' r " = {( u \r D ,j9Mr N ) I div( 7 V^) = in Q, u e H\n), 

supplier?) c T D }. 
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In the Calderon problem with partial data, the objective is to determine 
the conductivity 7 from the knowledge of C^ d,Tn for given (possibly 
very small) sets Td,Tn- Here d v is the normal derivative, and the 
conormal derivative ^fd v u\oo, is interpreted in the weak sense as an 
element of H~ l l 2 (dVt). Thus we have 

c v D ,r N c #i/2 (rD ) x H -^ 2 {Y N ). 

It is very useful to consider the related inverse boundary value prob- 
lem for the Schrodinger equation with partial data. There, the objec- 
tive is to determine a potential q G L°°{Vl) from the partial Cauchy 
data set 

C r q D ^ = {(u\r D , d v u\r N ) ; (-A + q) = in SI, u G H A (Q), 

Here, the space H A (Q) is defined by 

H A (Q) = {uE L 2 (VL) ; An G L 2 (Vl)}. 

It is known (see |BU01j ) that for u in this space, the trace u\qq and nor- 
mal derivative d u u\gn are in H~ l l 2 (dVl) and H~ 3 ^ 2 (dSl), respectively. 
Therefore 

Cr°< T » C H- l / 2 {T D ) x H-^ 2 (T N ). 
We formulate the relevant partial data problems as follows: 

Question 1. LetY d,Y ^ be open subsets of dQ and let 71,7% G L°°(Q). 
If 

°7i ~ °72 ' 

is it true that 7! = 72 ? 

Question 2. LetYo^N be open subsets of dQ and letqi,q2 G L°°(ft). 
If 

is it true that q± = q2? 

In both problems above, we think of u\au as Dirichlet data (volt- 
age) prescribed only on Yd, and we assume that one can measure the 
Neumann data of the corresponding solution (outgoing current) on Yjy. 
The Cauchy data set is determined by the Dirichlet-to-Neumann map 
(DN map) A 7 : H l / 2 (dfl) H- l / 2 (dfi), defined by 

A 7 : u\gn 1 — V 'fdvulan where u G H 1 ^) solves div(7Vw) = in fi. 

The partial Cauchy data set is obtained from the graph of A 7 as 

C t , d ' Tn = {(/|r D ,A 7 /|rJ; / e # 1/2 (<9fi), Bupp(/) C Y D }. 
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Also Cq D ' TN is obtained by restricting the graph of the DN map A q : 
u\an | — > d u u\on, where (—A + q)u = in Q, provided that is not a 
Dirichlet eigenvalue of — A + q in Q. 

One can think of three subcases of the above problems: 

• Full data: T D = T N = dfl. 

• Local data: = IV = T, where T can be any nonempty open 
subset of dfl. 

• Data on disjoint sets: and Tn are disjoint open sets. 

The most complete results are of course available for the full data 
case. Moreover, virtually all known results involve some version of the 
method of complex geometrical optics (CGO) solutions. This method 
has its origin in the works [Fa66], |Ca80j . and major results for full data 
inverse boundary value problems based on the CGO method are [SU87j . 
[HTllj in dimensions n > 3 and |Na96]. [AP06j . |Bu08] in the case 
n = 2. In particular, it has been proved that the set C® n,dCl determines 
uniquely a conductivity 7 G C 1 (fi) if n > 3 and a conductivity 7 G 
L°°(fi) if n = 2. See the survey |Uh09] for further references. 

For the partial data question where the sets or Tn may not be the 
whole boundary, we mention here four main approaches. Formulated in 
terms of the Schrodinger problem, it is known that C^ D,rN determines 
q in Q in the following cases: 

(1) n > 3, the set is possibly very small, and Ttv is slightly 
larger than dQ\TD (Kenig, Sjostrand, and Uhlmann |KSU07j ) 

(2) n > 3 and Y D = T N = T, and dfl \ T is either part of a hyper- 
plane or part of a sphere (Isakov |Is07j ) 

(3) n = 2 and = = T, where T can be an arbitrary open 
subset of dQ (Imanuvilov, Uhlmann, and Yamamoto |IUY10j ) 

(4) n > 2, linearized partial data problem, Td = Tn = T where 
T can be an arbitrary open subset of dfl (Dos Santos, Kenig, 
Sjostrand, and Uhlmann |D KSjU09| ) 

Here approach (1) requires roughly complementary sets T D and T N , 
whereas approaches (2)-(4) deal with the local data problem. Ap- 
proaches (l)-(3) also give a partial data result for the conductivity 
equation with the same assumptions on the dimension and the sets 
To, Tat. In (4), the linearized partial data problem asks to show injec- 
tivity of the Frechet derivative of A q at q = instead of injectivity of 
the full map q 1— > A q , when restricted to the sets T D and T^r. 
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It is interesting that each of the four approaches is based on a ver- 
sion of CGO solutions, but still the approaches are distinct in the sense 
that none of the above results is contained in any of the others. Ap- 
proach (1) is based on Carleman estimates with boundary terms for 
limiting Carleman weights, whereas approach (2) is based on reflection 
arguments and the full data methods of [SU87j . Approach (3) uses lim- 
iting Carleman weights with critical points and stationary phase, but 
involves complex analysis and is therefore restricted to two dimensions. 
Approach (4) is based on analytic microlocal analysis but so far only 
works for the linearized case. 

In a recent work [KS12j . we unified the approaches (1) and (2) and 
extended both of them, giving new partial data results in dimensions 
n > 3 also on certain Riemannian manifolds. Below, we will explain the 
approaches (1) and (2) from the new perspective obtained from [KS12] . 
and we will also give detailed proofs of certain partial data results in 
[KS12] restricting to the special case of Euclidean domains. We hope 
that the present treatment will be more accessible to readers familiar 
with Euclidean analysis than the geometric paper |KS12j . Approaches 
(3) and (4) give rather final results in the local data problem for n = 2 
and for the linearized case. We refer to the recent survey [GT13j on 
two-dimensional partial data problems for more details on approach (3). 
Approach (4) will be discussed in Section [7J and some open problems 
are listed in Section |HJ 

We list here some further references for partial data results, first 
for the case n > 3. The Carleman estimate approach was initiated in 
[BUOlj and [KSU07J. Based on this approach, there are low regularity 
results |Kn06] . |Zhl2j . results for other scalar equations |DKSjU07] , 
[KS07] [ChlT], [ CEE] an d systems [STTU] . stability results |HW06] . 
[HW07| . [Ti08] . |CDR12j . and reco nstruc tion results jNSTO] . The re- 
flection approach was introduced in [Is07j . and has been employed for 
the Maxwell system |COS09j, [Call] . Partial data results for slab ge- 
ometries are given in [LiUlO] , |KLU12] . A result analogous to Theorem 
12.71 was independently obtained in [IY12cj . 

In two dimensions, the main partial data result is |IUY10j which 
has been extended to more general equations [IUYllaj . combinations 
of measurements on disjoint sets [IUYllbj . less regular coefficients 
[IY12aj . and some systems |IY12b] . An earlier result is in [ALP05] . 
In the case of Riemann surfaces with boundary, corresponding partial 
data results are given in jCTllaj . jGTllbj . jAGTUllj . 

For piecewise analytic conductivities, uniqueness in the local data 
problem follows from boundary determination results |KV84] , |KV85] . 
See |Ge08| for a related construction. An early result for the linearized 
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problem in an annular domain in R 2 with no measurements on the in- 
ner boundary is in [Ha98]. In the case when the conductivity is known 
near the boundary, the partial data problem can be reduced to the full 



data problem [188], |AU04j . [Fl07j . [Be09] . |AK12j . |HPS12j . Also, we 



remark that in the corresponding problem for the wave equation, it has 
been known for a long time (see |KKL01j ) that measuring the Dirichlet 
and Neumann data of waves on an arbitrary open subset of the bound- 
ary is sufficient to determine the coefficients uniquely up to natural 
gauge transforms. Recent partial results for the case where Dirichlet 
and Neumann data are measured on disjoint sets are in |LO10j . |LQ12j . 

The paper is structured as follows. Section [TJ is the introduction. In 
Section [2] we give precise statements for the various partial data results 
in the literature, concentrating on the Schrodinger problem. Section [3] 
describes the main strategy for proving most of these results, by reduc- 
ing to a density statement for products of solutions that vanish on parts 
of the boundary. A Carleman estimate with boundary terms relevant 
for partial data results is proved in Section |H This is a special case of 
the corresponding estimate in |KS12j . but we give self-contained proofs 
in the Euclidean case. Section discusses the construction of CGO so- 
lutions vanishing on parts of the boundary, and Section E] explains the 
corresponding uniqueness results relying on the injectivity of a mixed 
Fourier transform/local attenuated geodesic ray transform of the un- 
known coefficient. The linearized problem is considered in Section [7J 
and Section [8] lists open questions related to the partial data problem. 

Acknowledgements. C.K. is partly supported by NSF, and M.S. is 
supported in part by the Academy of Finland and an ERC Starting 
Grant. 



In this section we give precise statements of the partial data results 
mentioned in the introduction. We will only consider the Schrodinger 
problem, since the conductivity problem can be reduced to that case 
by the well known relation 



This reduction works if the conductivities are sufficiently regular (one 
may additionally need to assume that the two conductivities agree on 
part of the boundary). 

It will be useful to introduce some notation. In all results below, 
we assume that O C 1" is a bounded open connected set with C°° 



2. Partial data results 



Kf = 7 



/ 2 A 7 ( 7 - 1 / 2 /) + 27^(^7)/, q, 



A( 7 1/2 ) 



1/2 
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boundary. If ip e C°°(Q) is a real valued function, we can write the 
boundary of Q as the disjoint union 

on = dtt+ U dVL- U dtt tSLn 

where 

dtt± = dQ±(cp) = {xedtt; ±Vip(x) ■ u(x) > 0}, 
<9fi tan = (9fi tan (v?) = {x E dn ; Vip(x) ■ v(x) = 0}. 

Here v is the outer unit normal of dQ. 

The functions ip that can be used in partial data results are typically 
limiting Carleman weights. This is a special class of functions, intro- 
duced in [KSU07j . that coincides with the set of harmonic functions 
(with some restriction on their critical points) if n = 2, and if n > 3 it 
consists of the following six functions up to translation and scaling: 

a ■ x, log 7— r^-, arg((a + i/3) • x), 



\x 



2 

2 



arg(e ty (x + iO), log 



\x-C\ : 



Here a and are orthogonal unit vectors, 6 € [0,2-71"), £ G M n \ {0}, 
and the argument function is defined by 

lm(z) 

arg(z) = 2 arctan — — . . , z<eC\(— 00, 0. 

v ' \z\ + Re(z) 

We refer to |DKSaU09] for more information and a thorough analysis 
of limiting Carleman weights that have no critical points. 

It is suggested by the methods of [KS12] that for a fixed limiting 
Carleman weight <p, measuring Neumann data on dQ + (ip) for Dirichlet 
data input of (</?), with no measurements required on dQ tan (<p), 
might be sufficient for determining the unknown coefficients (and this 
should also hold with dfl + ((p) and dQ-(ip) interchanged). All results 
described below can be understood in light of this idea, but most of 
them require measurements on dQ±((p) for a large family of different 
y?'s instead of just one (p. 



2.1. The result of [KSU07J. This result is stated in terms of the front 
and back faces of dQ with respect to some point xq which is outside 
the convex hull ch(f2) of Q. Note that if Q is strictly convex, the front 
face can be made arbitrarily small by placing xq close to the boundary, 
but in this case the back face will be very large. 



CALDERON PROBLEM WITH PARTIAL DATA 



7 



Theorem 2.1. Letfl cR n ,n> 3. If x e M n \ch(0) ; de/ine the front 
and back face of dQ by 

F(x ) = {i6ffl;(i- x ) ■ u(x) < 0}, 

B(xo) = {ieffl; (x — x ) • ^(x) > 0}. 

Lei r D ,rAr be open subsets of dfl with F(x ) C T D and B(x ) C T N . 
UQuQ2 e and if 

t/ien qi = q 2 - 

We also state another partial data result, which follows from the 
previous theorem by placing xq sufficiently far from Q. This is close to 
the earlier result of [BUOlj . 

Theorem 2.2. Let C M. n , n > 3. If a EM. n is a unit vector, define 
the front and back face of dfl by 

F(a) = {xedVt-a- v{x) < 0}, 

B(a) = {x edQ; a- u(x) > 0}. 

Let TdjTn be open subsets of dfl with F(a) cTj) and B(a) C T^. If 
qi,q2 € L°°(Q) and if 

/^r D ,rjv _ riT D ,T N 

then qi = q 2 . 

Both results can be understood as follows: partial measurements are 
sufficient for determining the potential provided that and IV satisfy 

dn^) u dn t M c r D , dn + &) u on t Uv) c 

for a suitable family of limiting Carleman weights. In the first theorem 
<p(x) = log | I where x' ranges over a small neighborhood of xq, 

and in the second theorem (p(x) = a' ■ x where a' ranges over a small 
neighborhood of a in the unit sphere. 

2.2. The result of |Is07| . This is a local data result where one 
measures both Dirichlet and Neumann data on the same set r = Tjj — 
Tn, but where the inaccessible part To = dfl \ F has to be part of a 
hyperplane. 

Theorem 2.3. Let fl C lR n ; n > 3. Assume that tt C {x n > 0}, let T 

be a nonempty open subset of dfl, and assume that T = dfl\T satisfies 
r C {x n = 0}. If qi ,q 2 e L°°{Q) and if 

r»r,r _ ^r,r 

then qi — q 2 . 
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By applying the Kelvin transform K(x) = x/\x\ 2 , this theorem im- 
plies a similar result (also proved in |Is07j ) where the inaccessible part 
of the boundary has to be part of a sphere. 

Theorem 2.4. Let Q C M. n , n > 3. Assume that Q is a strict subset 
of some ball B, let T be a nonempty open subset of dQ, and assume 
that T = dfl \ T satisfies T C dB. If q u q 2 G and if 

then qi = q 2 . 

These results can be understood so that measurements on T are 
sufficient to determine the potential if the inaccessible part Tq satisfies 
certain geometric conditions related to limiting Carleman weights. In 
the first theorem, it holds that 

T C p| <90tan(<A*) 
aeR n ,|a[=l 

a-e n — 

where <p a (x) = a ■ x and e n is the nth coordinate vector. Taking 
complements, this condition means that T should be sufficiently large 
in the sense that 

U (ffl + (^)u^y)cr. 

aeR n ,|a[=l 

a-en— 

In the second theorem, if the coordinates are normalized in such a way 
that B = i?(|e n , |), one has instead 

r c p| dn tan (K*<p a ) 

aeM",|a|=l 

a-e n — 

where K is the Kelvin transform and K*(p a is the limiting Carleman 
weight K*ip a (x) = a ■ x/\x\ 2 . 

2.3. The result of |IUY10] . This is a general local data result that 
is valid for two-dimensional domains. (It was extended to potentials 
with W X)V , p > 2, regularity in [IY12a J.) 

Theorem 2.5. Let (lei 2 be a bounded open set with C°° boundary, 
and let T be a nonempty open subset of dfl. If qi,q 2 £ C 4 ' a (Q) for 
some a > and if 

r»r,r _ ^r,r 

then q\ = q 2 . 
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The result is related to limiting Carleman weights as follows. Let T 
be any open subset of dfl, and let T = dQ \ T be the inaccessible part 
of the boundary. The proof of |IUY10] begins by showing that one can 
find a family of harmonic functions (limiting Carleman weights in 2D) 
{ip p } p( zs, where S is a dense subset of fl, such that 

T C p| dtt tan (Lp p ) 
pes 

and each ip p is a Morse function (its critical points are nondegenerate) 
having a critical point at p. We refer to |IUY10j and the survey [GT13] 
for more details about partial data problems in two dimensions. 

2.4. The result of |DKSjU09| . As it is explained in Lemma O 
below, if T D = T N = T, where T is a nonempty open subset of dfl, and 
if C T q f = C?f, then 

(qi - q 2 )u 1 u 2 dx = 

for any Uj G H&(fl) with (—A + qj)uj = in fl and supp^lan) C T. 
In |DKSjU09] , the authors consider a linearization of this assumption 
at qi = q 2 = 0. The main result of |DKSjU09l is: 

Theorem 2.6. Let Q be a bounded connected open subset ofM. n , n > 2, 
with connected C°° boundary. The set of products of harmonic func- 
tions in C°°(Q) which vanish on a closed proper subset F C dfl is dense 
inL\Q). 

While the proof of this theorem is quite different from the ones of 
the other results mentioned here, it does depend on the linear limiting 
Carleman weights <p(x) = a -x and the associated complex geometrical 
optics solutions e^', where ( 6 C" satisfies ( ■ ( = 0, or equivalently 
( = r(a + ij3) where a,^ e 1™ are unit vectors satisfying a ■ = 0. 
See Section [7] for a sketch of the proof of this result. 

2.5. The results of |KS12j . The partial data results in |KS12] are 
valid on a class of Riemannian manifolds of dimension n > 3, but here 
we will only mention some consequences for domains in M 3 . In all of 
the results below H C R 3 is a bounded open connected set with C°° 
boundary, and we consider the decomposition 

on = dn + (<p) u dn_(ip) u 9fi tan (^) 

with respect to a fixed limiting Carleman weight (p. We will also assume 
that dfltan is decomposed as 

dQ tan = T a U Ti 
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where T a is an open subset of dQt&n that is accessible to boundary 
measurements, and I\ is the inaccessible part. Further, we will assume 
that Tjj and are nonempty open sets in dfl with dQ_ UT a C Tp 
and dQ + UT a cT N . 

The content of the results below is that if the inaccessible part Tj 
satisfies a geometric condition, then Dirichlet measurements on dQ- U 
T a and Neumann measurements on dQ + UT a are sufficient to determine 
the unknown coefficient locally away from IY 

The first theorem corresponds to the case (f(x) = X\ and the case 
where the inaccessible part Tj is part of a cylindrical set. This result 
generalizes Theorem 12.21 (if one chooses E = 8Qq and Tj = 0) and 
Theorem 12.31 (if one chooses Qq and E so that fl fl {£3 = 0} C R x 
(dfl \E)). 

Theorem 2.7. Suppose that fl C R x flo where flo is a bounded open 
set with C°° boundary in R 2 . Let <^(x) = x\, and let E be an open 
subset of dflo such that Ti satisfies 

r, cMx (dtt \E). 

U G C{Vi) and if 

then q 1 — q 2 in Q fl (R x O) where O is the intersection ofQ with the 
union of all lines in R 2 that have <9f2 \E on one side. 

Note that if Qq is convex and E is connected, then the set O in the 
previous theorem is just the convex hull of E in Q and one recovers 
the potential in fl fl (R x ch^ o (E)). The next theorem is related to the 
logarithmic weight and generalizes Theorem 12.11 (choose E = dflo and 
Ti = 0). The inaccessible part of the boundary must now be part of a 
conical set. 

Theorem 2.8. Let Q C {£3 > 0} and let (p(x) = log|x|. Consider 
the hemisphere S 2 = {u G S 2 ; 0J3 > 0}, and let M be a compact 
subdomain of S 12 with C°° boundary. Let E be an open subset of dM 
such that Ti satisfies 

r i c{ru;r>0,u6 dM \ E}. 

If 1iil2 £ C{VL) and if 

then qi = q 2 in ft fl {ruj ; r > 0,u G 0} where O is the union of all 
great circle segments in such that dM \ E is on one side of the 
hyperplane containing the great circle segment. 
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The final theorem involves the weight <p(x) = arg(a;i +1x2) and does 
not have a counterpart in the previous results. It corresponds to a 
case where the inaccessible part of the boundary is part of a surface of 
revolution. 

Theorem 2.9. Let Q C {x± > 0} and let (p(x) = arg(xi +1x2)- Let 
S = {{xi, 0, X3) ; X\ > 0} be a half-plane with hyperbolic geodesies given 
by the half-circles (with R > and a G M.) 

(xi(t), X2(t), xs(t)) = (Rsint, 0, Rcost + a), t G (0, 7r) 

or the lines 

(x 1 (t),x 2 (t),x 3 (t)) = (t,0,a), t>0. 

Let Mo be a compact subdomain of S with smooth boundary, and let E 
be an open subset of dM such that Tj satisfies 

T i c{Re(dM \E);9e(-7r,7r)} 

where Rgx = (Rg(xi,X2),x 3 ) and Rq rotates vectors in M 2 by angle 9 
counterclockwise. That is, we assume that the inaccessible part Tj is 
contained in a surface of revolution obtained by rotating the boundary 
curve 8Mq \ E about the x^-axis. 
If 0.x 1 0.2 ^ C(Q) and if 

then qi = q 2 in Q fl {Rg(0) ; 9 G (— 7T, 7r)} where O is the union of all 
geodesies in S that have 8Mq \E on one side. 

Note that the above results are local results that allow to determine 
the unknown potential locally away from the inaccessible part IY The 
paper |KS12] also contains results where one obtains information on 
the potential near the inaccessible part Tj, but this information comes 
in the form of integrals along broken geodesic rays in the transversal 
manifold. In the Euclidean case, a continuous curve 7 : [0, L] — > Q is 
called a nontangential broken ray if 7 is obtained by following straight 
lines that are reflected in the standard way (angle of incidence = an- 
gle of reflection) whenever they hit 8Qq, all reflections are nontangen- 
tial, and all reflection points are distinct. One also needs a somewhat 
stronger assumption on the sets T D and T N (this assumption was made 
for simplicity and could be removed in many cases). 

Let us state a counterpart of Theorem 12 . 71 which involves the broken 
ray transform with exponential attenuation. Note that in the absence 
of the stronger assumption on To and T^, this result would generalize 
Theorem 12.71 
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Theorem 2.10. Suppose that f2 C R x Qq where fio a bounded open 
set with C°° boundary in R 2 . Let <^( x ) = x i> an d ^ E be an open 
subset of 8Qq such that l\ satisfies 



Let qi,q% E C(0), let T be a neighborhood of dQ + U dfl^ U T a , and 
assume that 



Then for any nontangential broken ray 7 : [0, L] — » Qq with endpoints 
on E, and given any real number \, one has 



Jo 

Here {■)" is the Fourier transform in the X\ variable, and qi — q 2 is 
extended by zero to R 3 \ Q. 

It follows from this result that one could recover the unknown po- 
tential also near Tj if one knew how to invert the broken ray trans- 
form, that is, to recover a function in Qq from its (attenuated) integrals 
over nontangential broken rays with endpoints on a given open subset 
E C 8Qq. This is a question of independent interest and there seem to 
be only partial results in this direction, see |Es04] . |I112j . 



In this section we give an outline of the proof of the partial data 
results described above. The proofs proceed in three steps: 

1. Via an integral identity, the partial data uniqueness question 
is reduced to showing that products of solutions of Schrodinger 
equations, which vanish on suitable parts of the boundary, are 
dense in 

2. Construction of a family of special complex geometrical optics 
solutions to Schrodinger equations, and showing that certain 
transforms of any function that is L 2 -orthogonal to products of 
complex geometrical optics solutions must vanish. 

3. Showing that the transforms arising in Step 2 are injective, 
which completes the proof that products of solutions are dense 
and also the uniqueness proof. 

This general outline roughly applies to all of the results in Section [21 
but the particulars (the choice of complex geometrical optics solutions, 
and the transforms in the final step) vary from case to case. 



r.cix (dn \E). 



wr,r _ ^r,r 




3. Strategy of proof 
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3.1. Reduction to the density of products of solutions. This 
step is achieved by the following basic integral identity. 

Lemma 3.1. IfT D ,T N C dQ are open and if C^^ Tn = C^ 2 d ' Fn , then 



/ (<?i - q2)uiu 2 dx = 
Jn 



for any Uj G H&(tt) satisfying (—A + qj)uj = in Q and 

suppfwil,^) C T D , supp(w 2 |an) C T N . 

Proof. Let Uj be as stated. Since C Fd ' Fn = C FD,rN , there is a function 
u 2 G H&(Q) with (—A + q 2 )u 2 = in Q, supp(w 2 |,9!ra) C T D , and 

(ui\r D ,d v ui\r N ) = {u 2 \r D ,d u u 2 \r N )- 
Using that u±, u 2 and u 2 are solutions, we have 

/ (<?i - Q2)uiu 2 dx = [(Aui)u 2 - ui(Au 2 )] dx 
Jn Jn 

= / [(A(ui - u 2 ))u 2 - {ui — u 2 )(Au 2 )) dx. 
Jn 

Now U\ —u 2 \on = 0, so in fact ui — u 2 G H 2 (VL) by the properties of the 
space H&(Q), see [BUOlj . Recall also that C°°(Q) is dense in H&(Q) 
and that u 2 \ dn G #" 1/2 (<9fi) and d u u 2 \ dn G #- 3/2 (<9ft). These facts 
make it possible to integrate by parts, and we obtain that 



(91 — q 2 )u\u 2 dx — / [{d v {ui - u 2 ))u 2 - {u\ - u 2 )(d u u 2 )} dS 
! Jon 

in the weak sense. The last expression vanishes since d v {ui— u 2 )\r N = 
and supp(M 2 |an) C T N . □ 

The next statement makes precise the reduction of the partial data 
problem to density of products. 

Corollary 3.2. Let VL dW 1 be a bounded open set with C°° boundary, 
let qi,q 2 G L°°(Q), and let and Tn be nonempty open subsets of 
d£l. If the set 

{u\U 2 ; Uj G H&(Q), (—A + qj)uj — in ft, 

suppfuil^) C T D , supp(u 2 |an) C T N } 

is dense in L 1 (f2), then the condition C^' Vn = Cg 2 D,rN implies q\ = q 2 . 



3.2. Complex geometrical optics solutions. To show the density 
claim in Corollary 13.21 one constructs special solutions Uj G H&(Q) of 



14 



CARLOS KENIG AND MIKKO SALO 



the equations (—A + qj)uj = in Q with the required boundary con- 
ditions. This will involve complex geometrical optics solutions having 
the form 

u x = e~ Tp (m 1 + n), 
u 2 = e Tip {m 2 + r 2 ) 

where r > is a large parameter, ip is a real valued weight function, 
mi and m 2 are complex amplitudes, and r± and r 2 are correction terms 
that are small when r is large. The solutions U\ and u 2 are chosen with 
opposite signs for the exponential [e~ TLp and e T(fi ) to make sure that 
these exponentials go away in the product, allowing one to consider 
the limit of U\U 2 as r — > oo. 

The point is to choose the amplitudes m 1; m 2 so that A(e _Tip m 1 ) rs 
and A(e T{p m 2 ) « in a suitable sense, ensuring that the approximate 
solutions e~ Tlf> mi and e T<p m 2 can be corrected to exact solutions of 
(—A + qj)uj = 0. The complex amplitudes often involve phase fac- 
tors e tT ^\ so the solutions contain complex exponentials e r<_</3+1 ^ 1 ) and 
e T(ip+v4> 2 ) _ This explains why these are called complex geometrical optics 
solutions. 

The choice of the weight function ip is crucial in this process. The 
correction terms r x and r 2 are obtained by solving the equations 

e T ^(-A + q 1 )(e-^r 1 ) = f 1 in Q, 

e -^(-A + q 2 )(e^r 2 ) = f 2 in ft 

where 

/i = -e rv (-A + g 1 )(e-^m 1 ), 

/ 2 = -e- T n-A + g 2 )(e^m 2 ). 

Thus one would like to have solvability results for the conjugated oper- 
ators e TLp {— A + q 1 )e~ Tip and e~ Tlp (— A + q 2 )e Tip that come with estimates 
such as 

lkj||i^(n) < C^WfjWL^n), 
for a constant C uniform over all sufficiently large r. Additionally, one 
needs some control of the boundary values of r± and r 2 on parts of the 
boundary. 

The above procedure can be carried out if (p is a limiting Carleman 
weight. We refer to [KSU07] . |DKSaU09] for more information about 
the microlocal condition characterizing these weights that ensures the 
right kind of solvability results. However, as discussed in Section [21 
limiting Carleman weights are just harmonic functions (with some re- 
striction on their critical points) if n = 2, and in R™ with n > 3 there 
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are only six possibilities up to translation and scaling. Moreover, as 
will be explained in Section H] below, one can expect relatively good 
control of the boundary values of r\ on dQ + (ip) and weak control on 
dQ tan (ip), and similarly relatively good control of the boundary values 
of r 2 on dQ-(ip) and weak control on <9fi tan (<y2) (this is due to the sign 
in the exponential e TT<p ). This will follow from a Carleman estimate (a 
weighted L 2 estimate with exponential weights depending on the large 
parameter r) with boundary terms. 

It remains to describe the construction of the amplitudes m\ and 
m 2 . The following three conditions are typical ones that one might like 
the amplitudes to satisfy: 

1. \\e^{-A+q 1 ){e-^m 1 )\\ L 2 (n) = || e -^(-A + g 2 )(e^m 2 ) \\ L 2 (n) = 
0(1) as r — )• oo, 

2- ||"^i|| J L2(an + ( ¥ >)) = \\ m 2 1| L 2 (an_( v )) = 0(1), and ||mj||i2( a n tail ( ¥ j)) = 
o(l) as r — )• oo, 

3. the set of limits {lim T _ 5 . 00 mim 2 } for all such mi and m 2 is a 
dense set in L 1 (fi). 

The first condition means that e~ Ttp m\ and e Tip m2 are approximate 
solutions in a suitable sense, allowing to find correction terms Tj with 
|l r j||L 2 (Q) = as r — > °° • The second condition comes from the fact 
that one wants roughly that 

uilan+uacita.n = 0> U2\dn.uan^n — 0; 

or equivalent ly 

r i|<9Q+u<9n ta n = ~ m i\dn + udn tan ^ r 2|an_uan t an = — m 2|an_U(9Qtan- 

Since the Carleman estimate with boundary terms will give relatively 
good control of the boundary values of rj on dQ + or dQ- and weak 
control on dQ t a,n, it is enough that rrij are bounded on dfl + or <9fi_ but 
they need to be small as r — > oo (or vanish completely) on <9f2 tan . The 
third condition follows since 

lim U\U2 = lim mim 2 . 

Thus, to show density of the products Uiu 2 , it is enough to show density 
of the set {linx^oo mxm^\ for all admissible mi and m 2 . This leads to 
the transforms considered in the final step. 

3.3. Injectivity of transforms. In the density of products approach, 
one needs to show that any function / in, say, C(Q), for which 



I fuiu 2 dx = 
n 
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for all solutions Uj G H&(Q) of (— A+qj)uj = in fl, must satisfy / = 0. 
By the discussion in this section, one can try to do this by constructing 
complex geometrical optics solutions with amplitudes rrij — rrij{- ; r, a) 
depending on the large parameter r and some additional parameters 
described by a which ranges over some parameter set A. Assume that 
for each a, there is a bounded measure \i a on Q such that 

lim m 1 (- ;t, a)m 2 { • ; r, a) = (j, a 

T— >O0 

in the weak topology of measures in Q. The orthogonality condition 
then implies that 

/ dfi a = 0, a e A. 
That is, the transform Tf of / vanishes, where 

Tf(a) = / fdfi aj a G A. 
Jn 

If this transform is injective, then / = and the density of products 
follows. 

Many different transforms have appeared in this connection. In the 
full data result of |SU87] . the transform T was just the usual Fourier 
transform. In the partial data result [BU.Qlj one obtained the Fourier 
transform in a cone, and also Isakov's approach |Is07] results in the 
Fourier transform. The partial data result |KSU07] involves a more 
complicated transform whose injectivity was proved by analytic mi- 
crolocal analysis, and this argument was simplified in |DKSjU07] where 
matters reduce to inverting the two-plane Radon transform. In the 
two-dimensional case, the transform appearing [Bu08j and |IUY10] is 
related to stationary phase. Finally, in |DKSaU09] , |KS12] a mixed 
transform involving the Fourier transform in one direction and X-ray 
transforms (or more generally attenuated geodesic ray transforms) in 
other directions was employed. 

4. Carleman estimates 

As discussed in Section El the first step in the construction of complex 
geometrical optics solutions with controlled boundary behaviour is a 
Carleman estimate with boundary terms. Here we give the proof of 
such an estimate for the conjugated operator e I/,//fe (— A + q)e~ tp / h in Q, 
where the limiting Carleman weight if is a linear function and h > 
is a small parameter (in the previous notation, h = 1/t). Following 
|KSU07j . it is useful to consider a slightly modified weight 

if £ = if + hf e 
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where f £ is a smooth real valued function in Q depending on a small 
parameter e but independent of h. The convexity of f £ will lead to 
improved lower bounds in terms of the L 2 (Q) norms of u and Vm. On 
the other hand, the sign of d„(p E in the boundary term of the Carle- 
man estimate will allow to control functions on different parts of the 
boundary. Of special interest is the set <9f2 tan where d u ip = 0, so one 
has 

d u ip e = hd u f e on <9fi t an- 

We would like to have d u f £ < on dQ tSL11 . It is not easy to find a 
global convex function f £ satisfying the last condition for a general set 
<9f2 tan . However, splitting f £ to a convex part whose normal derivative 
vanishes on dQ tan and another part which ensures the correct sign on 
dfltun will give the required result. We will use the notation D = —iV, 
(u,v) = (u,v)l^(q), \\u\\ = \\u\\l2^, and also L 2 inner products and 
norms on parts of dQ are indicated by a subscript. 

Proposition 4.1. Let Q C M n , n > 2, be a bounded open set with C°° 
boundary, let (p(x) = a ■ x where a G W 1 is a unit vector, and let k be 
a smooth real valued function in Q so that d u n = — 1 on <9f2. Let also 
q G L°°(Q). There are constants e,C ,h > with h < e < 1 such 
that for the weight 

h ^ 2 u 
(fi e = ip-\ — + flK 

e I 

where < h < h , one has 

h 3 h 2 

—(\d u (p £ \d v u, d„u)Qn_ {ips) + -^(Hl 2 + \\hDu\\ 2 ) 

< \\e^ h ((hD) 2 + h 2 q){e-^ h u)\\ 2 + h 3 (\d u ip £ \d u u, d v u) m+{ „ s) 
for any u G C°°(f2) with u\dn = 0. 

Proof. It is sufficient to consider the case a — e 1 , so (p(x) — x 1 and 

hx\ , 

ip £ = xi H — + tin. 

e I 

Let = e^/ h (hD) 2 e-^/ h , and write 

P Ve = A + %B 
where A and B are the formally self-adjoint operators 

A=(hD) 2 -\Vpe\ 2 , 

B = Vif £ ohD + hDo Vif £ = 2Vy? e • hD + -Ay? £ . 

i 
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Then, using the condition u\gn = 0, 

\\P^u\\ 2 = ((A + iB)u, (A + iB)u) 

= || Au\\ 2 + II^BwH 2 + i(Bu, Au) - i(Au, Bu) 

= || Au\\ 2 + \\Bu\\ 2 + (i[A, B]u, u) - ih 2 (Bu, d„u) dn 

= \\Au\\ 2 + \\Bu\\ 2 + (i[A, B]u, u) - 2h 3 ((d^ £ )d u u, d u u) dn . 

We compute the commutator 

i[A, B]u = h[{{hDf - |V^ £ | 2 )(2V^ e • Vu + (A<p £ )u) 
- (2V Ve ■ V + A<p £ )((hD) 2 u - \V Ve \ 2 u)} 

n 

= h [2V(hD) 2 ip £ ■ Vu + 4hD d k ^ £ ■ hDd k u + {{hD) 2 A^ £ )u 

k=l 

+ 2hDAip £ ■ hDu + 2V(f £ • V(\V(f £ \ 2 )u] 

n 

= h[4(ip"V<p £ ■ Vip £ )u - Ah 2 ^2 d jk if £ d jk u - 4h 2 VAip £ ■ Vu 

j,k=i 

-h 2 (A 2 <p £ )u]. 

Here f £ is the Hessian matrix of ip £ . Integrating by parts once, using 
again that u\qq = 0, yields 

(i[A, B]u, u) = 4h 3 (ip"Vu, Vu) + Ah{{y"Vip £ ■ Vy? £ )w, u) 

-h 3 ((A 2 pe)u,u). 

Now we use that (p is a linear function, so that tp" = |ei ® e x + hn" . 
Combining this with the previous computation for ||P^ £ u|| 2 , we obtain 



Ah 4 

\\P^u\\ 2 = \\M\ 2 + \\ Bu \\ 2 + — ll^i^ll 2 + 4^ 4 (k"Vm, Vu) 
Ah 2 

+ — \\(di<Pe)u\\ 2 + Ah 2 ((K"V(p £ ■ V Ve )u,u) - h 4 ((A 2 K)u,u) 

£ 

- 2h 3 ((d u (p £ )d u u, d u u) dn . 
We have V<p £ = V(f + h(^ £ - +Vk). Assume that h is so small that 

' (pVip 



h 



+ Vk 



< 1/2 on Q when < h < h . 
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Then 1/2 < |Vy? e | < 3/2 on Q, and for some constant C K > 1 

Ah 2 

\\P V M\ 2 > \\Auf + \\Bu\\ 2 + — \\hD lU \\ 2 - C K h 2 \\hDu\\ 2 

£ 

h 2 

^ IMI 2 ~ C K h 2 \\u\\ 2 - C K h 4 \\u\\ 2 - 2h 3 ((d v ip £ )d v u, d v u)do.- 

e 

We obtain a further lower bound for ||/iDm|| 2 from the term ||Au|| 2 : 
given any M > 0, and using that u\on = 0, one has 

h 2 \\hDu\\ 2 = h 2 ((hD) 2 u,u) = h 2 (Au,u) + h 2 (\Vif £ \ 2 u,u) 

From this lower bound for ||Au|| 2 and the trivial estimate ||-Bm|| 2 > 0, 
it follows that 

Ah 2 

\\P^u\? > 2Mh 2 \\hDu\\ 2 - M 2 h 4 \\u\\ 2 - 8Mh 2 \\u\\ 2 + ^L\\hD lU \\ 2 

- C K h 2 \\hDu\\ 2 + — \\u\\ 2 - C R h 2 \\u\\ 2 - C K h 4 \\u\\ 2 

s 

- 2h 3 ((d u ip £ )d l/ u, d u u) dn . 

Choosing M = (lGs)^ 1 , choosing e sufficiently small depending on C K 
and then choosing ho sufficiently small depending on e implies that 

\\P V M\ 2 > ^(\\u\\ 2 + \\hDu\\ 2 )-2h 3 ((d l/ (p £ )d„u,d„u) dn 

when < h < h Q . By further decreasing h we may add a potential 
q G L°°(Q) to get the estimate 

2||(P„ e + h 2 q)u\\ 2 > y^(IM| 2 + \\hDu\\ 2 ) - 2h 3 ((d u <p £ )d u u, d u u) dn 

Finally, we replace u by e x ^ 2e+K u, where u G C°°(Vt) and u\qq = 0, and 
use the fact that 

1/C < e x ' /2e+K <C onO. 
The required estimate follows. □ 

We now pass from ip £ to if in the boundary terms of the previous 
result, making use of the special properties of tp e on dQ. Note that 
the factor h 4 in the boundary term on {x G dQ ; —5 < d l/ ip(x) < h/3} 
below is weaker than the factor h 3 in the other boundary terms. This 
follows from the fact that d u ip £ = hd v K = —h in the set where d u ip 
vanishes, so one only has weak control near dQ tmi ((p). 
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Proposition 4.2. Let q G L°°(Q), and let <p{x) = ±xi. There exist 
constants C , h > suc/i that whenever < h < h and 5 > 0, one 
has 

o~h 3 h 4 h 2 

wit.. no 1 1 2 / 1 1 1 1 2 1 1 . 1 1 2 \ 

7^ 1 1 ^ u 1 1 -5}) +77" ll^ n llL 2 ({-<5<d^<V3}) + 7HIMI +||/ivw|| ) 

C-0 Oo Oo 

< \\e^\-h 2 A + ^ g )(e-^u)|| 2 + h s \\d u u\\h {{dvV > h/ 3 } ) 
for any u G C°°(Q) with u\gn = 0. 
Proof. Note that 

<9^ £ = ( 1 + -(p) d v ip + hd v K — ( 1 + -<£> ] ct<£> - h. 



We choose ho so small that whenever h < ho, one has for x G 

1 , h , N 3 

On the set where d v tp(x) < —5, we have 

|0„¥> e | > 5/2. 
If —5 < d u ip < h/3, we use the estimate 

\d^ £ \ > h/2. 

Moreover, \d u ip £ \ < Co on dQ. Since {d v <p < h/3} C {d u (p e < 0} and 
[d v p e > 0} C {c^ > /i/3}, the result follows from Proposition 14.11 
after replacing Co by some larger constant. □ 

We can now obtain a solvability result from the previous Carleman 
estimate in a standard way by duality (see [BU011 IKSU071 INSlOj ). 
There is a slight technical complication since the solution will be in L 2 
but not in H l . To remedy this we will work with the space 

H A (tt) = {u G L 2 (VL) ; Au G L 2 (tt)} 

with norm ||w||# A = |M|l 2 + ||^ m ||l 2 - It is known (see |BU01j ) that 
H A (Q) is a Hilbert space having C°°(fi) as a dense subset, and there is 
a well defined bounded trace operator from H A (Q) to H^^^dVL) and a 
normal derivative operator from H A (Q) to H^/ 2 {dVt). We also recall 
that if u G H A (Q) and u\qq G # 3 / 2 (<9fi), then u G H 2 {Vt). 

Proposition 4.3. Let Q C M. n be a bounded open set with C°° bound- 
ary, let q G L°°(fi) ; and let ip(x) = ±X\. There exist constants 
Cq, t > such that when r > r and 5 > 0, then for any f G L 2 (Q) 
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and /_ G L 2 (S_ U So) there exists u G L 2 (Q) satisfying e T ^u G H&{VL) 
and e Tip u\QQ G L 2 (<9f2) such that 

e~^(-A + q)(e^u) = f into, e^u\ s _ uSo = e^/_, 

and 

i/ere S± and So are the following subsets of dVt: 

S- = {d vV < -5}, S = {-5 < d uV < l/(3r)}, S + = {d vV > l/(3r)}. 

Proo/. Write Ln = e T ^(-A + q)(e~ T ^v) and r = 1/h, r = l//i - We 
rewrite the Carleman estimate of Proposition 14.21 as 

(<5r) 1/2 ||^n|| L 2 (5 _) + \\d u v\\ L 2 {S()) + r\\v || + ||Vn|| 

^Coll^H+Cor 1 / 2 !!^!!^^). 

This is valid for any 5 > 0, provided that r > r and u G C°°(f2) with 
^|an = 0. 

Consider the following subspace of L 2 (Q) x L 2 (S + ): 

X = {(Lv, d v v\ s+ ) ; v G C°°(Ti), v\ dn = 0}. 

Any element of X is uniquely represented as (Lv, d„v\s + ) where v\qq = 
by the Carleman estimate. Define a linear functional I : X — > C by 

l(Lv,d u v\ s+ ) = (v,/)l3(o) - (duV,f-) L *(S-L>s )- 
By the Carleman estimate, we have 

|Z(Lu,3„v|s + )| < ||v|| ll/H + II^Hl^s.)!!/-!!^^) 

+ ||Mz,2(s )||/-||l2(5o) 

< ^(r-l/H + (5r)- 1 / 2 || / _|| i2(s _ ) + ||/_|| i2(So) ) 
x (ULull+r^lia^ll^)). 

The Hahn-Banach theorem implies that I extends to a continuous linear 
functional T : L 2 (tt) x t' 1/2 L 2 (S+) C such that 

Hill < CbCr-l/ll + (5r)- 1 / 2 ||/_|| L2(5 _ ) + HZ-IU^)). 

By the Riesz representation theorem, there exist functions u G L 2 (Q) 
and u + G L 2 (S+) satisfying 1(10,10+) = (w,u) L 2 {n) + (w + ,u + ) L 2 {s+) . 
Moreover, 

IMU a (n) + t~ 1,2 \\u + \\ L 2 {s+) 

< Co(r- l \\f\\ + (Sr)-^\\f4ms^ + Wf-W^so))- 



22 CARLOS KENIG AND MIKKO SALO 

If v G C°°(fi) and v\qq = 0, we have 

(4.1) (Lv, u) L 2 (n) + (d u v, u + ) L 2 (s+) = (v, f) L 2 {Q) 

- (d u vJ_) L 2 (s _ uSo) . 

Choosing v compactly supported in Q, it follows that L*u = f, or 

e~ TV {-/\ + q){e TV u) = f in Q. 

Furthermore, e Tip u G H^(Q). 

If w,v G C°°(Q) with = 0, an integration by parts gives 

(Lv,w) = -{e~ Tip d v v,e TLp w) L 2 {dQ) + (v,L*w). 

Given our solution u, we choose Uj G C°°(Q) so that e Tlfi Uj — > e^u in 
H&(p,). Applying the above formula with w = «j and taking the limit, 
we see that 

(Lv,u) = -(e~ T(p d v v, e T(p u) L 2(pa) + (v,L*u) 

for v G C°°(Q) with v\gn = 0. Combining this with ( 14. ip . using that 
L*u = /, gives 

{.d u vJJ) L 2 {s _yj So) + (d u v,u + ) L 2 {s+) = (e~ Tlp d u v,e Tlp u) L 2 {m) . 
Since c^u can be chosen arbitrarily, it follows that 

e r</, w|s_us = e Tip f-, e Tip u\ s+ = e Tlp u+. 
We also see that e TLp u\ dn G L 2 (dVt). □ 

5. Complex geometrical optics 

We will now describe a construction of CGO solutions that is relevant 
for the proof of Theorem 12.71 Suppose that Q C R x Q where f2 is a 
bounded open set with C°° boundary in M 2 . Let <p(x) = xi, and let E 
be an open subset of dQ such that Tj satisfies 

r.cix (dn \E). 

Let also g G L°°(fi). We wish to construction a solution w G H&(Q) of 
the equation 

(-A + q)u = in tt 

where 

w = e -^i( m _|_ rQ ) 
and s is a slightly complex large frequency, 



s = t + i\ 
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where r > (eventually r — > oo) and A 6 M is fixed. The use of a 
slightly complex frequency instead of a real frequency allows to intro- 
duce another real parameter A in the CGO solutions, which makes it 
possible to employ the Fourier transform in the x\ variable. 
Inserting the ansatz for u in the equation, we need to solve 

where 

/ = - e SXl {-A + q){e- SXl m) = -(-A + 2sd x - s 2 + q)m. 

It will be useful to look for an amplitude m independent of x\, so 
m = m(x') where x = (x\,x') and x' = (x2,Xs) G M. 2 . Then / has the 
simpler form 

/ = (-A x > - s 2 + q)m. 
We seek for an amplitude m e C 2 (Q ) satisfying 

\\(-A x , - s 2 )m\\ L 2 (Uo) = 0(1), ||m|| L2(no) = 0(1), 
and the boundary values should satisfy 

IMU 2 (i?) = 0(1), \\m\\ L 2 {mo \ E) = o(l) 

as t — > oo. These conditions have been chosen to be compatible with 
Proposition 14.31 and they can be interpreted so that m should be an 
approximate eigenf unction, or quasimode, of the Laplacian in fl with 
frequency s. In fact, we will arrange so that 

m\dn \E = 0. 

If we can find such an m then Proposition 14.31 together with the fact 
that the inaccessible part Tj satisfies Tj C M x (8Qq \ E) will allow to 
find a correction term r with ||r||£2(n) = o(l) as r — > oo. 

Let us now describe one construction of such an m. We choose a 
straight line 7 in M. 2 that intersects f2 but does not meet <9f2 \ E, and 
will construct an amplitude m that concentrates on 7. We look for m 
in the form 

m( x >) = e is ^ x,) a(x'). 
A computation gives that 

(-A x ,-s 2 )m = e isqp [s 2 (\V?p\ 2 - l)a - is(2Vip ■ Va + (A x ,ip)a) - A x ,a] . 

It is enough to choose ip and a so that the following equations are valid 
in VLq\ 

|Wf = l, 2V^ • Va + (A x ,ip)a = 0. 
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The first equation is an eikonal equation, and distance functions are 
solutions. Choose a ball Cl with Q CC Cl , and choose some point 
x' Q G Qq \ that lies the line 7. Define 

if){x') = \x' — x' \, x' G fi . 

Then if) G C°°(fio) satisfies \Vip\ = 1- The second equation above is a 
transport equation for a, and has the solution 

a(r,6) = r - 1/2 b(6) 

where (r, 9) are polar coordinates in R 2 with center at x' Q , and b is any 
function in C°°(S' 1 ). (Note that if>(r, 9) = r in these coordinates.) The 
point is that if the line 7 is given in the (r, 9) coordinates by r h-> (r, 9q), 
then choosing b supported very close to 9 and independent of r will 
result in the boundary conditions 

IM|i2( B ) = HdfMi? = 

as r — > 00. 

Combining the above amplitude construction with Proposition 14.31 
results in the following existence result for CGO solutions. (See |KS12t 
Section 6] for the full details of the proof.) 

Proposition 5.1. Suppose that Q C R x Q where Q is a bounded 
open set with C°° boundary in R 2 . Let <^(x) = ±x± and decompose 
dfl tail (<f) as T a U Ti where the closed set r, satisfies, for some open 
subset E of OQq, 

r.clx (dtt \E). 

Let also q G L°°(Q). Choose a ball flo with Qq CC Clo, let x' G f2o\f2ci; 
and let (r, 9) be polar coordinates in R 2 with center at x' Q . Let 9q G S 1 
be such that the line r \-t (r, 9 ) meets £l but not dilo \ E. There exists 
a solution u G H&(Q) of the equation 

(-A + q)u = inVL 

having the form 

u = e- sv {e isr r-^ 2 b{9) + r ) 
where b G C°°(S' 1 ) is supported very close to 9q, 

\\ro\\iP(a) = o(l) 

as t — > 00, and 

supp(u\ dn ) C r 
ifT is any open set in dQ with dfl T ((f) U T a C T. 
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6. Uniqueness results 



We now describe how to complete the outline given in Section [3] and 
prove Theorem 12.71 using an injectivity result for a certain transform 
(in the end of the section we discuss briefly the proofs of Theorems 
I2.8H2.10)) . If / is a piecewise continuous compactly supported function 
in R 3 , the relevant transform of / is given by 



where A G R and 7 is a line in R 2 . That is, we are taking the Fourier 
transform of / in the X\ variable with frequency 2A and then taking 
the attenuated X-ray transform, with constant attenuation — 2A, along 
lines in R 2 ,. We choose the parametrization j(t) = au±_ + too, where 
u> G S 1 is the direction vector of the line 7, u± is the counterclockwise 
rotation of u by 90 degrees, and a is the signed distance to the origin. 
(The choice of parametrization will not be too important below.) 

Proposition 6.1. Suppose that Q C R x Qq where Qq is a bounded 
open set with C°° boundary in R 2 . Let <p(x) = x\, and let E be an open 
subset of dQo such that Tj satisfies 



r.cMx (dfl \E). 
Ifqi, q 2 G C(n), ifdQ_ UT a cr D and dQ + U T a C T N , and if 



for all A G R and for any line 7 in R 2 that does not meet dQo \ E. 
(Here, we extend q± — q 2 by zero to .) 

Proof. By Lemma 13.11 our assumptions imply that 



for any Uj G Ha(Q) satisfying (—A + 51)1*1 = (—A + q 2 )u2 = in Q 
and 

supp(wi|an) C T D , supp(u 2 |an) C T N . 

Fix A G R, let s = r + iX where r > is sufficiently large, and fix a 
line 7 in R 2 that does not meet dfl \ E. We use Proposition 15.11 to 




91 <?2 



then 



T(gi-g 2 )(A, 7 ) = 




2(» 
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find solutions u\ and u 2 satisfying the above conditions and having the 
form 

Ul = e- SXl (e isr r-^ 2 b(e)+n), 
u 2 = e SXl (e isr r- 1/2 b(6)+r 2 ) 



n 

OO POO 



where (r, 9) are polar normal coordinates in R 2 whose center is outside 
of VLq and such that the line 7 is given by r ^ (r, 9 ), b G C7°°(S' 1 ) is 
independent of r and supported near 9 , and 

\\rj\\ = o(l) as r — > 00. 

We then have 
= lim / (qi — q 2 )uiu 2 dx 

(qi ~ q2)(x u r,9)e- 2lXxi e- 2Xr \b{9)\ 2 d9 dr dx x 

-00 Jo 

Here we have extended qi — q 2 by zero outside Q. The result follows by 
choosing b to approximate a delta function at 9q. □ 

The uniqueness result now follows from the Helgason support theorem 
for the X-ray transform [He99], which states that if / is a piecewise 
continuous compactly supported function in R 2 that integrates to zero 
over any line that does not meet a compact convex set K C M 2 , then 
fW 2 \K — 0. The idea of reducing the attenuated X-ray transform to 
the unattenuated one comes from |DKLS13] . 

Proof of Theorem \2. 7[ Denote by / the extension of q\ — q 2 by zero to 
1R 2 , and let 

00 



00 

2 



Let also K be the convex hull of dVt® \ E in R . By Proposition 16. 1[ it 
is enough to show that if 

T/(A >7 ) = 

for all A G R and for any line 7 that does not meet K, then f(x±, x') = 
whenever x% G R and x' G R 2 \ K. 

The condition Tf(X, 7) = implies that 



e- 2M f(2X, 1 (t))dt = 

for any line 7 that does not meet K. Setting A = 0, the Helgason 
support theorem implies that f(0,x') = for x 1 G R 2 \ K. We now 
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differentiate the above identity with respect to A and set A = to 
obtain 



;-2t)/(0, 7 (t)) + 2^(0, 7 (t)) 



dt = 0. 



But we already saw that f(0,x') = for x' G M 2 \ K. The Helgason 
support theorem then gives that 

df 

-f (0,x) = for x G M 2 \K. 
Repeating this argument shows that 

r\ \ k 

' f(0, x') = for a;' G M 2 \ AT 



for all > 0. But x\ h-> f(xi,x') is compactly supported, hence £i i-> 
/(^i, x') is analytic, and we obtain that x') = for all d 6l and 
x' G 1R 2 \ K. The result follows upon inverting the Fourier transform 
in x\. □ 

Theorems 12.81 and 12.91 can be proved by the same general argument 
as above. However, the fact that one has nonlinear limiting Carleman 
weights (the log weight and arg weight) leads, after a suitable conformal 
scaling, to a situation where the original domain Q with Euclidean met- 
ric is replaced by a compact Riemannian manifold (M, g) with smooth 
boundary, and the condition 

n c r x n 

is replaced by the condition 

(M,g)c(RxM ,g), g = c(e®g Q ) 

where c is a positive function (conformal factor), and the transversal 
domain Q with Euclidean metric is replaced by a transversal manifold 
(Mo, go) that is compact with smooth boundary. Also the integrals over 
lines in fl are replaced by integrals over geodesies in (Mo, go). This 
setup is similar to the results for the anisotropic Calderon problem in 
[DKSaU09j . |DKLS13j . For the log weight, (M ,g ) will be a domain 
in the hemisphere, and for the arg weight (Mo, go) will be a domain 
in hyperbolic space. Correspondingly, one uses geodesies in the sphere 
and in hyperbolic space. In both cases, if (Mo, go) is chosen to have 
strictly convex boundary, then it is a simple manifold with real-analytic 
metric and analogues of the Helgason support theorem are available. 
The full proofs are given in [KS12J. 
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The proof of Theorem 12.101 requires a more general complex geomet- 
rical optics construction than the one described above. The reason is 
that one wants to get information on integrals along broken rays that 
may touch the inaccessible part Tj, and one then needs amplitudes in 
the complex geometrical optics solutions that concentrate along these 
broken rays and are very small on the inaccessible part. Such ampli- 
tudes may obtained from a Gaussian beam quasimode construction. 
This construction in connection with complex geometrical optics so- 
lutions is employed in [DKLS13] for non-reflected geodesies (full data 
case) and in |KS12j for reflected geodesies (partial data case), and we 
refer to these papers for the details. 

7. The linearized case 

Here we sketch the proof of Theorem 12.61 A first reduction is to show 
that it suffices to prove a local version of the theorem (see |DKSjU09 



Section 2]). The global version follows from the local one by using 
ideas in the spirit of the Runge approximation theorem, developed in 
an unpublished work of Alessandrini, Isozaki, and Uhlmann. Thus, 
matters reduce to proving the following "local" version. 

Theorem 7.1. Let Q C M. n , n > 2, be a bounded connected open set 
with C°° boundary. Let xq £ dQ and let F be the complement of some 
open boundary neighborhood of x . There exists 5 > such that given 
any f £ L°°(f2) ; if we have the cancellation property 

/ fu\U2 dx = 
Jn 

for all Uj £ C°°(f2) with Auj = in Q and Uj\p = 0, then f = in 
B(x Q ,8) n Q. 

The next step is to use conformal transformations (in particular 
Kelvin transforms) of harmonic functions to reduce to the following 
situation: x = 0, Q C {x £ M. n ; \x + e%\ < 1} where e\ = (1, 0, . . . , 0) 
is the first coordinate vector, the tangent hyperplane to is given 
by {xi = 0}, and F C {x £ dQ ; x\ < —2c} for some c > 0. See 
DK SjU09[ Section 3] for this reduction. 



From this point on, the proof is inspired by the proof of Kashiwara's 
"watermelon theorem" in analytic microlocal analysis. We introduce 
the Segal-Bargmann transform of a function / £ L°°(IR n ) with compact 
support by 

Tf(z)= [ e-^ z ^ 2 f(y)dy, z = x + i£eC n . 
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The Segal-Bargmann transform is related to the microlocal analysis 
of analytic singularities of a distribution. We mention the a priori 
exponential bound 

(7.1) |r/(^)|<(27r/ i r/ 2 e^l Im ^||/|| L o 

and the fact that if / is supported in {x\ < 0}, we can improve this to 

(7.2) \Tf(z)\ < (27r/i)" /2 e^ |Im2|2 - |Re ^ 12 



L°° 

when Re^i > 0. Both of these bounds are straightforward. Note also 
that when we multiply by (2n)~ n ^ 2 , when z = x G M™, we obtain the 
Gaussian heat kernel, and hence if / has compact support we have 

(7.3) lim {27rh)- n/2 Tf{x) = f{x) for a.e. x G R n . 

h— >0 

The strategy is to use the cancellation property (we extend / to be 
outside of fi) to obtain an exponential decay of Tf when x G M. n 
is close to 0, thus yielding the vanishing of / near 0, as desired. In 
order to accomplish this, let us temporarily consider the case n = 2 
and define 7 = ie% + e 2 = (i, 1) G C 2 . Note that {7, 7} is a basis of C 2 , 
and that {CgC 2 = is the union of two complex lines C7 U C7. 
Hence, it is easy to see if e > is small enough, then any z G C 2 with 
\z — 2ie\\ < 2e may be decomposed as a sum of the form 

(7.4) z = ( + r], (.( = r] -r 1 = 0, 

where \( — j\ < Ce, \r) + j\ < Ce. This last fact extends to C n , n > 2, 



see |DKSjU09j Section 3]. 



Recall that the exponentials e~ l ~ , £ G C n ( ■ ( = 0, are harmonic 
functions. We need to modify these by adding correction terms to 
obtain harmonic functions u satisfying the boundary requirement u\f — 
0. Let x £ C^°(R n ) be a cutoff function which is supported in {x\ < 
—c} and equals 1 on {x% < —2c}. Consider the solution w — w( ■ ; Q 
of the Dirichlet problem 

■ Ox 

Aw = in Q, w\ 9Q = -e 1 h x \au- 
We have the following bound on w: 

(7-5) IklU^n) < c i\\ e ~ l ~x\\m/z(dn) 



where Im^i > 0, from the choice of \ and our normalization of Q. 
We have the cancellation property 

f(x)u(x,C)u(x,r])dx = 0, (■( = i 1 -r] = 0, 



n 
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where u(x, Q = e~ % ~ +w(x, (), which is a harmonic function in C°°(il) 
and satisfies u\p = 0. The identity above can be rewritten as 



f(x)e~i xit:+v) dx 



- f f(x)e 



f(x)e hX ^w(x,7])dx 



hX ' r 'w(x, () dx — / f(x)w(x,()w(x,ri)dx. 



Using the estimate (17. 5p . we obtain 



f{x)e~^ +v) dx 



< c 3 \\f\\ LOO (i +h- i \r,\y/ 2 (i+h- i \c\)- l/2 

X min { Im Cl >Im ?7! } e ^ ( | Im C' I + 1 Im ??' | ) 



when lm(i > 0, Im^ > 0, and ( ■ ( = rj ■ 7] = 0. In particular, if 
|C — a l\ < Cea, \r] + 07] < Cea, e < 1/(2C), then 



f(x)e-^ x < c+v) dx 



< dh- 1 



L°°(f2) 



ca ICea 

e 2h e h . 



Hence if z G C™ and \z — 2iae±\ < 2ea, e small enough, using a rescaled 
version of the decomposition 17.41 gives 



(7.6) 



f(x)e h x ' z dx 



< C A hT x 



L°°(f2) 



ca 2 C z a 

e h e h . 



To relate the last estimate to the given estimates on the Segal- 
Bargmann transform of /, we use the well-known formula 

e -^(z-y) 2 = e -i(2irh)- n/2 J e-£e-* y < t+lz) dt 

which gives 

Tf(z) = (2nh)- n / 2 




e 2^ 



(z 2 +t 2 ) e -iy(t+iz) 



f(y) dtdy. 



For our /, supported in Q and verifying the cancellation of the integral 
in our hypothesis, the estimate (17. 6p and the formula above allow us 
to improve the estimate (17. 2p . If Re^i > 0, then 



\Tf(z)\ < (2nh) 



-n/2 



e ^(|Im^| 2 -|Re 2 | 2 -t 2 ) 



< 



e ^-(|Im 2 | 2 -|Re Z | 2 ) 
(27Th) n / 2 



e 2h 



. J\t\<m 



e-rM t+ ^f(y)dy 



dt 
dt 



e 2h 



t\>ea 



e h 



hv(t+iz) 



f(y) d y 



dt 
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Using (17. 6p with z replaced by t+iz, when — 2aei | < ea and \t\ < ea 
we obtain 



(7.7) \Tf(z)\ < e 27r(l Im ^l 2 -l R ^l 2 ) SU p f 

It I <ea J 



-^ t+x) f(y) dy 



+ V2e^ ez %-^ J\f(y)\dy 



< C.hT 1 



7±-(|Im 2 | 2 -|Re;z| 2 ) 
L™(Q)e- 2h 



ac 2Csa £ a ea 

e 2h e h + e ih e h 



provided that \z — 1ae\\ < ea. Choosing e < c/(8C) and a > (c + 
4e)/(e 2 ), we obtain the bound 

\Tf(z)\ < 2C 5 h- 1 \\f\\ L ~ {n) e^ I ™ z ^ Rc ^-f). 
Combining (17. ip . (17. 2p and ( 17. 6p we have 



where x' G 



on— 1 



and 



when Re Z\ < 0, 



(Imzi) 2 - (Re^i) 2 , when Re^i > 0. 



Now, inspired by the proof of the "watermelon theorem", we use the 
following: 

Lemma 7.2. Let b,L > 0. Let F be an entire function in C such that 

e fc F(s) < <^ , r| , 

[ e 2t, ty/iere | s — L\ < b. 

Then for all r > 0, there exist c',5 > (independent of F) such that 

\F{s)\<e-^ + ^ 

when | Res | < 5 and \Jms\ < r. 

The proof of the lemma rests on a harmonic majorization argument, 
which exploits the subharmonicity of — (Im s) 2 + (Re s) 2 (see DKSjU09 



Lemma 4.1 and Remark 4.2] for the details of the proof). Next we apply 
Lemma 17.21 to 

hTf(s,x>) 



F{s 

We obtain in particular that 



C||/IU°°(n) 

\Tf(x)\ < Ch-'WfU^e-^, xeQ, \ Xl \ < 5 
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for 5 small. Multiplying by (2rrh)~ n ^ 2 and letting h — > 0, using (17. 3p . 
we deduce that f(x) = for x G f2, — 5 < x\ < 0. The completes the 
proof of Theorem 17.11 and hence of Theorem 12.61 

8. Open problems 

The following is a list of some questions related to the partial data 
problem that are open, as far as we know, at the time of writing this. 
The first question concerns the local data problem for n > 3. (As 
discussed above, this result for n = 2 is known at least for sufficiently 
regular coefficients.) 

Question 1. (Local data in dimensions n > 3) If Q is a bounded 
domain in IR n , n > 3, if T is an arbitrary nonempty open subset of dQ, 
and if gi, g 2 G L°°(Q), show that C^f = C^f implies gi = g 2 . 

The next question concerns data on disjoint sets. The general case is 
open even for the linearized problem in any dimension. Partial results 
for n = 2 are given in }IUYllbj . 

Question 2. (Data on disjoint sets in dimensions n > 2) If 

is a bounded domain in M™, n > 2, if and are arbitrary disjoint 
open subsets of dQ, and if q x , g 2 G L°°(£l), show that C^^ N = C^ 2 D > rN 
implies qi = q 2 - 

One can also ask for optimal regularity conditions for the coefficients 
in partial data results. If n = 2 the full data result is known for L°° 
conductivities |AP06j . 

Question 3. (Local data for nonsmooth conductivities) If f2 is 

a bounded domain in R 2 , if T is an arbitrary nonempty open subset of 
dtt, and if 71, 72 G L°°(fi) satisfy 71, 72 > c> 0, show that C^f = C£.r 
implies 71 = 72. 

The next question concerns optimal stability for partial data results. 
In the full data case it is known that in general one has a logarith- 
mic modulus of continuity for determining an unknown coefficient from 
boundary measurements, and this result is optimal. We refer to the sur- 
vey |A107j . Stability results for partial boundary measurements based 
on the reflection approach [HW07] , [Call] and in the case where the co- 
efficient is known near the boundary [AK12] also have logarithmic sta- 
bility. However, available results for the Carleman estimate approach 
[HW06| . |Tz08] . |CDR12] seem to involve weaker moduli of continuity 
(log log or worse), and one can ask if logarithmic stability still holds. 
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Question 4. (Optimal stability for the Carleman estimate ap- 
proach) What is the optimal stability for the partial data uniqueness 
result of |KSU07j ? 

Finally, there are many open questions related to partial data for 
elliptic systems. In the introduction we mentioned results for systems 
when n = 2 and for the reflection approach when n > 3, but it seems 
that there are no partial data results for systems using the Carleman 
estimate approach in n > 3 except for [ST10J. As an example, one can 
consider the time-harmonic Maxwell equations for which the full data 
result is known [OPS93J. 

Question 5. (Carleman estimate approach for systems) Prove 
an analogue of the partial data result of |KSU07j for the time-harmonic 
Maxwell equations with scalar coefficients as in [OPS93J. 
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